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1. Introduction 

Heavy flavour production at hadron colliders has been the subject of extensive the¬ 
oretical and experimental studies for more than twenty years. The discovery of the 
top quark has offered an excellent opportunity to test QCD predictions much more 
reliably than in the case of bottom or charm, thanks to the smaller value of a s 
and the relatively minor impact of long-distance effects, the top having no time to 
hadronize before decay. At present, all comparisons between theory and data concern 
tt pair production; a crucial role in the satisfactory agreement between predictions 
and experimental results is played by the next-to-leading order (NLO) QCD cor¬ 
rections [1, 2, 3, 4], which enlarge the leading-order cross section by about 30% at 
the Tevatron. A companion process to pair production is that in which a single top 
quark is present in the final state. In such a case, a weak-interaetion Wtb vertex is 
involved, and thus the single-f cross section is smaller than the one for tt (in spite 
of being favoured by phase space volume), which so far has prevented observation 
of such a production mechanism by Tevatron experiments. In terms of Standard 
Model physics, single-f production is a direct probe of the weak interactions of the 
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top, which in fact constitutes the main interest of single-f signals. Amongst other 
things, this may lead to measurements that have not been performed so far, namely 
of the CKM matrix element V t b, and of the b parton density. Single-f production is 
in addition an important background for many searches for new physics, and can in 
general be seen as an effective way to study new physics phenomena in the heavy 
sector. 

For single-f searches, or counting experiments in which single-f is a background, 
it is crucial to have a reliable estimate of the number of events expected, i.e. of the 
total rate. In this respect, NLO results are mandatory, also in view of the fact that 
they allow a sensible assessment of the size of unknown contributions of higher orders. 
Calculations of fully-differential NLO single-f cross sections have been performed in 
refs. [5, 6, 7, 8] and, including NLO top quark decay, in refs. [9, 10, 11, 12, 13]. On 
the other hand, in order to optimize acceptance cuts in an experimental analysis, or 
to perform full detector simulations, one needs realistic hadron-level events, which 
are obtained with Monte Carlo event generators that incorporate the simulation of 
parton showers and hadronization models. 

The complementary benefits of fixed-order computations and parton shower sim¬ 
ulations have been discussed at length in the literature, as well as the advantages of 
combining them into a framework which would retain the strong points of each of 
them. The MC@NLO approach [14, 15] (we shall refer to these papers as to I and 
II respectively hereafter) provides a way of achieving this, by allowing one to match 
cross sections computed at NLO in QCD with an event generator. No modifications 
to the latter are necessary, and therefore existing parton shower Monte Carlos can 
be used for this purpose. 

Although the MC@NLO formalism has been defined in full generality in I, ex¬ 
plicit implementation details have been given there only for processes with no final- 
state QCD emissions at the level of hard reactions. Such a case has been considered 
later in II, with the implementation of ff and of bb production. In the context of 
MC@NLO, a process-independent calculation is required for each type of soft and/or 
collinear singularity which appears in the NLO real matrix elements. A quick inspec¬ 
tion of the processes implemented so far (see ref. [16]) should convince the reader that 
the only singularity structure untreated is the final-state collinear one. We shall deal 
with this singularity in the present paper. It must be clear that, as for all of the other 
singularities which have been studied previously, our formulation will not depend on 
the fact that the specific single-f production process is considered here: in the deriva¬ 
tion of the analytical formulae the nature of the hard reaction is irrelevant (which is 
further evidenced by the fact that the inclusion of single-f production in MC@NLO 
relies significantly on results obtained in I and II). What we achieve here is there¬ 
fore, besides the addition of an important process to the MC@NLO framework, the 
capability of including other processes in MC@NLO without the need of performing 
further analytical computations, notably those having final-state (massless) partons 
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at lowest order. 

The paper is organized as follows: in sect. 2 we discuss single-f production in 
the context of fixed-order computations and Monte Carlo simulations. Sect. 2.1.1 
reviews the status of the matrix elements used in the present computation. We limit 
ourselves here to implementing the s- and f-channel production mechanisms, and 
neglect spin correlations in production. In sect. 2.1.2 we show that some changes 
can be made in the subtraction formalism [17, 18] upon which MC@NLO is based, 
which leave its analytical expression unaffected, but improve its numerical stabil¬ 
ity. We then proceed to sect. 2.2, where we write down the approximate singlc-t 
production cross sections generated by HERWIG, which enter the definition of the 
MC subtraction terms needed for the matching with NLO results. Implementation 
details of MC@NLO, concerning in particular the simultaneous presence of initial- 
and final-state collinear singularities, are given in sect. 3. We present results for 
singlc-t production at the Tevatron in sect. 4; phenomenological studies, including 
results for the LHC, will be the subject of a future paper. Finally, conclusions and 
future prospects are reported in sect. 5. Some technical details are collected in the 
Appendices. 

2. Single-top cross sections 

Each process in MC@NLO is based on two main building blocks: a fully-exclusive 
NLO computation; and the knowledge of the so-called MC subtraction terms, which 
are closely related to the first non-trivial order in the formal a s expansion of the 
HERWIG Monte Carlo result. We shall treat these two issues in turn. 

2.1 NLO computation 

Fully-exclusive observable predictions do not strictly exist in QCD: the theory has 
finite resolution power, in the sense defined by the KLN theorem. However, we can 
conventionally talk of fully-exclusive computations, as those in which the cancella¬ 
tion of the infrared singularities is formally achieved analytically in an observable- 
independent manner, and the four-momenta of all of the final-state partons are avail¬ 
able for defining the observables - this does not violate the KLN theorem, since the 
formal cancellation mentioned above actually occurs only in the case of infrared-safe 
observables. Fully-exclusive computations are crucial for the matching of NLO cross 
sections with Monte Carlos, since the latter need to know the four-momenta of all the 
particles involved in the hard process in order to compute the initial conditions and 
the various branching probabilities for the parton showers. Modern computations of 
this kind are based on universal subtraction or slicing formalisms; we shall discuss 
the one used within MC@NLO in sect. 2.1.2. Before doing that, we give some details 
specific to the matrix elements for single-f production. 
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2.1.1 Matrix elements 


The lowest-order parton level processes are customarily divided into three classes 
that will also serve to categorize the NLO contributions. They are shown in fig. 1. 
In the first diagram the single top quark is produced in the annihilation process 

u + d —■> f + b , ( 2 . 1 ) 

via a time-like W boson, and is therefore called the s-charmel process. In the second 
diagram, the initial bottom quark is converted into a top quark via the exchange of 
a kb-boson 

b + u —»f + d , ( 2 . 2 ) 

and is therefore called the f-channel process. The final two graphs represent the Wt 
process in which the top quark is produced in association with a real kb 

b + g —> kb + t. (2.3) 

The cross section for this process occurring at the Tevatron is very small and we 
neglect it in this paper. For the LHC this process becomes non-negligiblc however. 
Note that in reactions (2.1) and (2.2) we have only listed the CKM-dominant com¬ 
binations of quark flavours, but all CKM-allowed combinations are included in this 
paper. Consistently, the b quark is always assumed to be massless. 



Figure 1: Leading order diagrams for single-f production in the (1) s-channel, (2) f-channel 
and (3) kbf-mode. The f-quark line is doubled. 


In NLO one must include virtual and real corrections to the s- and f-channel 
processes. The virtual corrections consist of vertex corrections to diagrams (1) and 
(2) in fig. 1, together with the self-energy corrections to the f-quark line 1 . We shall 
not discuss these corrections in detail, nor give their explicit expressions, as these 
are already given in the literature. To prepare a remark on single-antitop production 
further below, we recall here that the vertex correction in the first diagram of fig. 2 
is proportional to the lowest order vertex y M (l — 75 ) because only light quark lines 
are attached to it. If the top quark line is attached as in the second diagram of fig. 2, 
a second form factor appears at NLO, proportional to the difference ($ — p^/rrit. 
A similar situation occurs in the f-channel. 

graphs vanish since they involve a single colour matrix on a fermion line, i.e. a null trace. 
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Figure 2: Virtual vertex corrections to s-channel single-t production. 


Concerning the real-emission corrections, we categorize these processes by the 
dominant CKM contributions, as follows 


ud —> tbg , (2.4) 

ub —-> tdg , (2.5) 

bd —■> tug , (2.6) 

ug —■> tbd , (2.7) 

dg —> tbu , (2.8) 

bg —■> tdu . (2.9) 


A rather detailed discussion of these processes, and how they are assigned to s- and 
f-channel, can be found in the next section. 





Figure 3: Diagrams contributing to ug —> tbd. 


The calculation of the single-7 cross section is perfectly similar to that for the 
singlc-f described above, after charge conjugation. It may be perhaps less apparent 
that the second vertex form factor, mentioned above, proportional to (pi — pb)/rn t 
remains unchanged, since the quark propagators change the sign of the mass term. 
However the charge flow of the conjugated amplitude is also reversed, resulting in an 
unchanged expression. 

2.1.2 Subtraction procedure 

In order to implement a process in MC@NLO, its NLO cross section must be com¬ 
puted according to the subtraction formalism presented in refs. [17, 18] (denoted as 
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(1) (2) (3) (4) 

Figure 4: Diagrams contributing to bg —> tdu. 

FKS henceforth). The basic idea in FKS is that of partitioning the phase space of 
the final-state partons involved in real-emission contributions, in such a way that 
the resulting regions do not overlap, cover the whole phase space, and each of them 
contains at most one collinear and one soft singularity. In each of these regions it 
is natural to select the one parton (called the FKS parton here) with which the sin¬ 
gularities are associated. Denoting by the generic real matrix elements, this 

amounts to writing 2 

i = Z) s i 0)+ Z!‘ s « ) ’ < 2 ' 10 ) 

i ij 

M (r) = ^sj 0) M (r) + • (2.11) 

i ij 

The FKS parton is labelled with i in eqs. (2.10) and (2.11). The first term on the 
r.h.s. of eq. (2.11) gives a divergent contribution (i.e., a contribution which has to 
be subtracted) only in the infrared regions in which parton i is soft and/or collinear 
to one of the initial-state partons. Analogously, the only infrared regions in which 
the second term on the r.h.s. of eq. (2.11) is divergent are those in which parton i is 
soft and/or collinear to final-state parton j. More precisely, denoting by p a and k a 
the four-momenta of the initial- and final-state particles respectively, we have 


= u < 212 ) 

lim S$ 0) = 1, (2.13) 

Mpi 

lim = 1, (2.14) 

lim4 1) = l, (2.15) 

ki\\kj 


2 The notation of refs. [17, 18] has been slightly changed here in order to simplify the discussion. 
Functions S of the present paper play the same role as functions 0 in ref. [18]. 
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while all the other infrared limits not explicitly listed above are zero 3 . Eqs. (2.12)- 
(2.15) are the only properties of the S functions used in the analytical computations 
of refs. [17, 18]; their actual functional forms away from the infrared limits are only 
relevant to numerical integrations. It should be stressed that all partons in the final 
state may induce a divergence of the real matrix elements; to take this fact into 
account, the role of FKS parton is given to each parton in turn, which is formally 
expressed in eq. (2.11) by the sum over i that appears on the r.h.s. there. 

After the phase space of the final-state partons is effectively partitioned through 
eq. (2.11) into different infrared-singular regions, FKS chooses a different phase- 
space parametrization in each of these regions. It must be clear that the phase space 
is always the same, i.e. that relevant to the n particles involved in real-emission 
processes; the only difference between the various regions is in the choice of the in¬ 
tegration variables which are left after getting rid of the 5 functions that appear in 
the phase-space definition. The integration variables are chosen to perform the nec¬ 
essary analytical integrations in an easy way, and to facilitate importance sampling 
in numerical integrations. The key variables in the phase-space parametrization as¬ 
sociated with 7y 0) are the energy of parton i (directly related to soft singularities), 
and the angle between parton % and one of the initial-state partons (directly related 
to initial-state collinear singularities). For the energy of parton i and the an¬ 
gle between parton i and parton j (related to a final-state collinear singularity) are 
chosen instead. Obviously, the indices % and j are dummy here (phase spaces are 
flavour blind), and therefore there are only two independent functional forms for 
phase spaces in FKS, which loosely speaking are relevant to initial- and to final-state 
emissions. More details, and specific functional forms, are given in appendix B. 

After the partition of the phase space, achieved by means of S,j 0) and it is 
the matrix elements that determine whether a singularity actually occurs in a given 
region of such a partition. As a general rule, one should choose the simplest possible 
forms for the S functions that still allow subtraction of all singularities. Although 
this is by no means mandatory (a region without singularities will simply give a 
finite contribution to the cross section), it is beneficial for well-behaved numerical 
computations. Since singlc-t matrix elements have a singularity structure much sim¬ 
pler than that of the matrix elements considered in refs. [17, 18], the S functions will 
also be simpler here. We also want to use the present process as a test case, and 
will define the S's as smooth functions of invariants, at variance with the original 
formulation of refs. [17, 18], in which they have been expressed as products of 0 
functions. 

We start by denoting the four-momenta entering an NLO tree-level single-t pro- 

3 The superscripts (0) and (1) are legacy notation from ref. [17], where these S'-functions are 
related to jet-finding algorithms, and the superscripts indicate the algorithm step at which a merging 
takes place. 
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duction process as follows 


a(pi) + Pip*) —* t(h) + 7 (^ 2 ) + S(k 3 ) , (2.16) 

where a and f3 are the incoming partons from the left (jp\ > 0 ) and from the right 
(p 2 < 0) respectively; 7 and 5 denote light final-state partons. We shall use the 
following shorthand notation 

(a,fct,p,5) (2.17) 

for the momentum assignment of eq. (2.16). 

We first consider process (2.4); the treatment of processes (2.5) and (2.6) is 
identical 4 . We assign momenta as follows: 

(u, d; t, b, g ). (2.18) 


By inspection of the relevant Feynman diagrams, we immediately conclude that the 
only singularities are associated with the gluon: the final-state light quark cannot give 
rise to a collinear divergence, being in all cases connected to a W boson. Therefore, 
for such processes the gluon will always be the FKS parton and, according to the 
discussion given at the beginning of this section, we can choose the S functions in 
such a way that the only non-zero ones are and ■ In particular, with the 
following forms 


,(o) = _ (h-ki) a (k 3 -k 2 ) a _ 

(k 3 -k 1 ) a (k 3 -k 2 ) a + (k 3 -pi) a (k 3 -p 2 ) a 

,(i) = _ (k 3 'Pi) a {k 3 - p 2 ) a _ 

' 32 (k 3 -k 1 ) a (k 3 -k 2 ) a + (k 3 - Pl ) a (k 3 -p 2 ) a 


(2.19) 

( 2 . 20 ) 


equations (2.10)-(2.15) are fulfilled (i = 3, the gluon being the FKS parton). In 
eqs. (2.19) and (2.20) a is an arbitrary positive real number; the physical results 
will not depend on a, and their stability against the variation of a will constitute a 
check of the correctness of our implementation. It is clear that the numerators of 
eqs. (2.19) and (2.20) will act as damping factors for final- and initial-state collinear 
singularities respectively; the larger a, the stronger the damping. Formally, in the 
a —» cx) limit we could recover the ©-based implementation of the S functions of 
refs. [17, 18]. More pragmatically, we shall use the freedom in the choice of a to 
improve, if necessary, the numerical stability of the result, and will study its impact 
on the number of negative-weight events in MC@NLO, see section 4. 

We now turn to the case of process (2.7); the corresponding Feynman diagrams 
are shown in fig. 3. There are only initial-state collinear singularities in this case, 
due to the splittings g —> bb (graph 2) and g —*■ dd (graph 3). On the other hand, 

4 It is immediate to see that the procedure adopted here to disentangle the singularities of 
(a, /?; t, 7 , S) works identically for (/?, a; t, 7 , (5). 



these two diagrams do not interfere: the former contributes to the t-channel cross 
section, the latter to the s-channel one. Since s- and f-channel contributions are 
integrated separately, we are in the same situation as process (2.4) (i.e., only one 
parton can give singularities), except for the fact that no final-state singularities are 
present in this case. Therefore, we can set = 0 here, which implies S (0> = 1. It 
also implies that we are in the same situation as that treated in I (which also applies 
to many other processes implemented in MC@NLO). This situation now naturally 
appears as a particular case of a more general implementation in which singularities 
are disentangled by means of S functions. 

Since process (2.8) is completely analogous to process (2.7), we finally deal with 
process (2.9), whose Feynman diagrams are shown in fig. 4. Of those, graphs 1 and 4 
contribute to the Wt mode, which has not been considered here and are therefore 
dropped, while graphs 2 and 3 contribute to the t-channel. Graph 2 (graph 3) is 
singular when the u ( d ) is emitted collinearly to the initial-state gluon; since the two 
diagrams do interfere, we disentangle the singularities by means of the S functions. 
We assign the momenta according to 

(b,g;t,d,u), ( 2 . 21 ) 


and we define 

s (o) = _ ( k r p 1 ) a {k r p 2 ) a _ 

2 {k 2 -pi) a (k 2 -p 2 ) a + ( k 3 -p 1 ) a (k 3 -p 2 ) a ’ 

S (P) = _ {k 2 -pi) a {k 2 -p 2 ) a _ 

3 ( k 2 ■ Pi) a ( k 2 • p 2 ) a + ( k r p 1 ) a (k r p 2 ) a ’ 


( 2 . 22 ) 

(2.23) 


which again fulfill equations (2.10)-(2.15). Although the same arbitrary parameter 
a as in eqs. (2.19) and (2.20) has been used here, this is in fact not necessary; we 
could introduce another free parameter, independent of a. 

We conclude this section by stressing that the functional dependences of the 
S functions given above are correlated with the momentum assignments chosen for 
the corresponding subprocesses. For example, eqs. (2.22) and (2.23) imply that the 
FKS parton will have four-momentum k 2 and k :i respectively. Clearly, the subtrac¬ 
tion formalism is independent of the particular labeling adopted for each process. 
Therefore, through a relabeling we can always assign four-momentum k :i to the FKS 
parton. Such relabeling is a purely formal trick to render manifest the local matching 
between NLO matrix elements and MC subtraction terms. 

As far as processes (2.4)-(2.6) are concerned, we pointed out before that only 
the gluon can play the role of FKS parton. Thus, the momentum assignment in 
eq. (2.18) and the analogous ones 


(u, b; t, d, g) (2.24) 

(b, d] t, u, g) (2.25) 
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are what we want; as a consequence, the S functions have still the forms given in 
eqs. (2.19) and (2.20). For the process in eq. (2.7), see fig. 3, we noted that in 
the cases of s- and f-channel contributions the singularities arise from the splittings 
g dd and g —> bb respectively. Therefore, we assign momenta as follows 


(■ u,g;t,b , d) 

s — channel, 

(2.26) 

(u, g\ t, d, b) 

t — channel. 

(2.27) 

and analogously for process (2.8) 

( d,g-,t,b,u ) 

s — channel, 

(2.28) 

(d,g;t,u,b) 

t, — channel. 

(2.29) 


Finally, owing to the fact that S^ «-»• when /c 2 <->• & 3 , we write 

M {r \b,g]t,d,u) = Si 0) M {r) {b,g-t,d,u)+Si 0) M {r \b,g-,t,d,u) 

— M^ r \b,g-,t,d, u) + M^ r \b, g-,t,u,d) . (2.30) 

In other words, we shall assign the momenta in process (2.9) in two different ways 

(■ b,g;t,d,u ), (2.31) 

(■ b,g]t,u,d ), (2.32) 

and for each of them we multiply the corresponding matrix element times given 
in eq. (2.22); as shown in eq. (2.30), this is fully equivalent to eqs. (2.21)-(2.23). 

As a final remark, we note that when keeping the same ordered notation (2.17) 
after charge conjugation the treatment of real emission corrections to anti-top pro¬ 
duction is perfectly analogous, and the inclusion of single-7 requires no extra work. 

2.2 MC cross sections expanded to NLO 

As discussed in I and II, in order to construct the MC subtraction terms one needs 
the cross section obtained by keeping the first non-trivial order in the a s expansion of 
the parton shower Monte Carlo that will be matched with the NLO computation. As 
in the previous papers, the explicit results presented here are relevant to HERWIG. 
The most general form of the MC cross sections is given in eq. (II. 5.1), which we 
rewrite as follows 

da da l L ' l) ’ ( 2 - 33 ) 

MC - - - MC 

fd L l 

where the index // generically indicates a collection of labels which unambiguously 
identify the 2 —> 3 partonic subprocess. The index L assumes the values +, —, / 1 , 
and /2 (the latter two were denoted by Q and Q in II). The index /, which differs 
per colour structure, assumes the values q t -q v where q % and q 3 are the four-momenta 
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of the colour partners relevant to the emission considered; in this way, the shower 
scale is 

K = |i| = la-ait (2.34) 

In II we had / = s,t,u (and Eq = |/|/ 2 ), but in the case of unequal masses this 
is not a convenient notation. Equations (II.5.2)-(II.5.5) are also unchanged apart 
from notation 


da^’ l) 

h 1 

dv { ~' l) 

dcrl fl ’ l) 

da (M 

H' 


MC 

MC 

MC 

MC 


-(iy fa Hl) ( X U/ z + ) fb H2) (®2i) daj+V 


,(0 


fa Hl \xii)fi H2 \x 2 i/z W ) dal 


f^\x lf )f^\x 2f )da^ 

f^\x lf )f^\x 2f )da^ 


dx 


1 / 


dx 


1 / 


dx^i dx 2 i , 

MC 


dxn dx 2 i , 

MC 

dx 2 f , 
dx 2f , 


(2.35) 

(2.36) 

(2.37) 

(2.38) 


where the flavours a and b of the incoming partons depend on the value of /i. The 
short-distance cross sections that appear on the r.h.s. of eqs. (2.35)-(2.38) can be 
read from eq. (II.5. 6 ) and eq. (II.5. 8 ) 


dd^ 


da^ 


«s d£±_ 
2vr ^ 

^djl 

27r tW 

^ J a 





(2.39) 

(2.40) 


where the @’s account for HERWIG dead regions (see sect. 4.3 of II), and the flavours 
a', b ', and the values of /i' can be determined by considering the possible collinear 
splittings of the corresponding NLO tree-level processes. 

As in II, we use unbarred and barred symbols to denote quantities relevant to 
2^3 and 2^2 processes respectively. The momentum assignments for the former 
are given in eq. (2.16), while for the latter we use 


a'(pi) +/3\p2) —> t(ki) + 7 \k 2 ), 


(2.41) 


which we shorten in a way similar to eq. (2.17) 


(«', t, i) 


(2.42) 


In MC cross sections expanded to NLO, 2^2 momenta (entering da on the r.h.s. 
of eqs. (2.39) and (2.40)) are obtained by means of a suitable projection of the 
corresponding 2 —> 3 momenta. The exact form of the projection is specific to 
the parton shower MC matched to the NLO computation, and for HERWIG can be 
worked out as was done in II. Here, we need to extend the formulae given in II, in 
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(u, d ; t, b ) 

(d, u) t, b) 

(n, dr, t, b, g) 

±(pi-P2);fi,2(ki-h) 


( d,u;t,b,g) 


±(pi -P2); fiaih-h) 

('u,g-,t,b , d) 

-{Pi 'Pi) 


( d,g-,t,b , u) 


~{Pl-p2) 

(g,u-,t,b, d) 


+ {pl-p2) 

C g,d',t,b,u ) 

+ {Pl ’Pi) 



Table 1: Short-distance contributions to MC subtraction terms, for the s-channel. The 
two columns correspond to the two possible .s-channel Born cross sections, distinguished by 
the direction of the incoming partons. For a given process, the entries show the emitting 
legs, and in round brackets the value of the shower scale Eq (up to a sign), according to 
the possible colour flows. 



( 6 , u) t, d) 

( 6 , d] t, u ) 

{b, u; t, d, g) 

+Ji(Pi-h); -J2{p2-h) 


(■ b,d-,t,u,g) 


+,fi(prhy, -,f2{p2-k 2 ) 

(b,g\t,d, u) 

~{p2-h) 


(. g,u-,t,d,b) 

+{pi-h) 


( b,g-,t,u , d) 


-{p2-h) 

(. g,d]t,u,b ) 


+{pi-k i) 


Table 2: As in table 1, for the t-channel, with bu- and M-initiated Born processes. 



(■ u , b ; t, d) 

(d, b ; t, u ) 

(u,b;t,d,g) 

+J2(pi-h); -Ji{p2-ki) 


{d, b; t, u, g) 


+,f2(prhy, - Ji{p2-k 1 ) 

(u,g-,t,d,b) 

-(p2-h) 


(g,b;t,d, u) 

+{pi-k 2 ) 


{d, g\ t, u, b) 


~{p2-ki) 

(g,b;t,u, d) 


+{pi-k 2 ) 


Table 3: As in table 1, for the f-channel, with ub- and d6-initiated Born processes. 

order to treat the case of final-state partons with unequal masses; explicit results are 
given in appendix A. 

As far as flavour combinations are concerned, it is simpler to read eqs. (2.39) 
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and (2.40) from right to left, since this follows the logic which forms the basis of 
a parton shower. The MC starts with a Born-level (2^2 for single-f production) 
process, and then lets each leg branch in all kinematically- and flavour-allowed con¬ 
figurations possible. This implies that several 2^3 processes may be generated 
starting from a given 2 —» 2 process. We list all such processes explicitly in tables 1- 
3; the non-void entries give non-zero contributions to eq. (2.33). Thus, the index /x 
that classifies the 2^3 partonic processes can simply be chosen so as to count all of 
the quantities that appear in the first columns of the tables. Parton legs where the 
branchings occur are denoted by +, —, f \, and f 2 (/i always coincides with the top 
quark); given the parton that branches, and the hard subprocess, a colour connection 
is established which fixes the shower scale Eq unambiguously. The shower scales to 
be used in eqs. (2.39) and (2.40) are equal to the absolute values of the dot products 
listed in tables 1-3. We finally point out that the momentum assignments for the 
2^3 processes in the tables above are the same as those adopted (after relabeling) 
in the context of the pure NLO computation. This puts the NLO and MC cross 
sections on the same footing from a notational viewpoint, which will be convenient 
for the formal manipulations to be carried out in the next section. 


3. MC@NLO 


The definition of the MC@NLO formalism is given in eq. (1.4.22) or eq. (II.2.1): 


X ~vi 


^ J dxidx2d<j)3^J r 1 


(3) 

MC 


dE 


(/) 


d<p 3 


dE, 


■jri 2 ) 

•J MC 


dE 


■(/) 


1 /dE 


d<t> 3 

<c+) 


d(p 2 dx 


dE„ 
-1- a 

ct d(j) 3 

(c-) 


+ 


d(j) 3 

1 /dE 


MC 


T, 


MC 

(b) 


+ 


dE 


(sv) 


+ 


dE 




d<p 2 dx 


1 /dE 


d(j) 2 dtp 2 

(c+) 


xX d(p 2 dx 


dE 


(c-) 


ct + d<p 2 dx 


Ct 


(3.1) 


There is a minor difference of notation with respect to eq. (II.2.1): the indices for the 
sum over all partonic processes are denoted in the present paper by n, consistently 
with what was done in sect. 2.2. We refer the reader to I and II for all details 
relevant to the formalism. Singlc-t production is the first process implemented in 
MC@NLO in which both and SjP are non-zero for certain i and j. This has 
direct implications for eq. (3.1), which we now discuss. 

As shown in sect. 2.2, for a given choice of the index /x which classifies the partonic 
processes the radiation pattern in the MC cross section is determined by the values 
of the indices L and l. On the other hand, the possible radiation patterns at the 
NLO level are determined by the S- 0) and functions. Inspection of sect. 2.1.2 
and of tables 1-3 shows that, for a given /x, there are at most one and one 
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<S !I) functions which are non-vanishing. Formally, this corresponds to defining two 
single-valued functions i(fj) and j{/i) such that and may be different 

from zero. This allows us to define the following quantities 


c(in) _ c(0) 


o(out) _ <j(l) 


(3.2) 


where the labels IN and OUT are to remind us that and SjP select kinematics 
configurations relevant to initial- and final-state collinear emissions respectively. Note 
that one of the functions in eq. (3.2) may be still be vanishing (which is the case 
for in processes (2.31) and (2.32)), but there cannot be other non-vanishing S 
functions. In any case, from eqs. (3.2) and (2.10) we obtain 


S<™> + 5< OUT) = 1 V/i. 


(3.3) 


Note that, since we have exploited relabeling invariance to assign the four-momentum 
k 3 always to the FKS parton, we have f(/i) = 3. Furthermore, since the other 
massless final-state parton has four-momentum /e 2 , we also have j (/i) = 2. However, 
eq. (3.2) holds independently of relabeling invariance. Furthermore, it is clear that 
an analogous equation must hold for any kind of hard reaction, and not only for 
singlc-f production, as a direct consequence of the definition of the FKS partition. 

The term is proportional to the real-emission matrix elements (see eq. (1.4.12) 
and eq. (1.4.13)); thus, according to eq. (2.11), in eq. (3.1) we understand 


>v</) 

d(j) 3 


,yd/) 

_ c(in) 

M d(p3 

.yK/) 

c(in) 


.yK/) 

I c(out) 

** #3 

rlT U) 

, c(out) I * 

“ d4°" T) 


(3.4) 

(3.5) 


where we have introduced two different parametrizations {d(f)^ } and d(f> 3 ° UT ' > ) of the 
three-body phase space d(p 3 , analogously to what is done in FKS in the context of 
pure NLO computations (see sect. 2.1.2). Their explicit forms, which are irrelevant 
in what follows, will be given in app. B. 

In eq. (3.1) each point (x\,X 2 , 4> 3 ) corresponds to a 2 —> 3 kinematic configuration 
(called HI). In previous MC@NLO implementations, a definite 2^2 configuration 
(called S) was chosen given (aq, x 2 , 4> 3 ), according to a mapping whose form is 

dictated by HERWIG. The definition of such an unique mapping requires elaborate 
manipulations of the MC subtraction terms since, as shown in eq. (II.B.32) and 
eq. (II.B.33), initial- and final-state emissions would naturally lead to the definitions 
of two different mappings and Following the same arguments as in 

app. B of II, we could implement single-f production using an unique Vm^s] as 
discussed there, however, this may degrade the numerical accuracy in the integration 
step and the unweighting efficiency. Furthermore, the two mappings and V^g 1 
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are a perfect match to the FKS phase-space partition, which enters eq. (3.1) through 
eq. (3.5); the mapping PEs (EE?) can be naturally expressed in terms of the 
variables of d^E (dcj) EE- In practice, we replace eq. (1.4.23) and eq. (1.4.24) with 
the following pairs of equations 


and 


4 N, = E/ 


ry(f) 

dxxdxzd^ f ST ] —hj 
d4 3 


WL=±) 


4 IN) = E 


dx\dx 2 d(j )E 


ev 

. .dsj? 

c(in) M 
J , (in) 


d4 iN) 


d(Pz 


+ 


dE 


MC 

<L=±) 


ct 


#E 


5, 


(!N) , dS ( 6 ) 


x 2 


d<j). 


lo 


(in) 
2 

,y( c +) 


+ 


i^( sv ) 


d(j). 


^dx 


d(j), 

+ 


+ 


(in) 

2 

,yr(c-) 


Ct 


d(j). 


^dx 


To 


Ct 


>v=r( c +) 


d 0 ! 




+ 


,wr( c -) 


d<p 


^dx 


j h OUT) = E [ dxidx 2 d(f)^ >VT ^ ( 5^ out) 


,y(/) 

^,7 


d 0 ( OUT) 


jv(i=/l,2) 


d(p 


(out) 


MC 


r(0UT) 


21 ^ 


dxidx 2 d(p y 3 


(out) 


_ p(out) 


d4 /j 


#e t) 


+ 


dE? /l,2) 


ct 


5, 


(out) 


lo 


dE 


(ft) 


#E T) 


+ 


dE 


M OUT) 

(«f) 


d4 OUT) 


in such a way that the total rate is now 


^ _ d IN ) I r( IN ) I /-(out) , A out) 

°tot - + % + J s + % 


(3.6) 


(3.7) 


(3,8) 


(3.9) 


(3.10) 


which replaces eq. (1.4.25). It should be clear that eqs. (3.6)-(3.9) are a direct 
consequence of the definition of MC@NLO: in fact, eq. (3.1) is recovered by insert¬ 
ing Ec on the r.h.s. of eqs. (3.6) and (3.8), and 4c on the r.h.s. of eqs. (3.7) 
and (3.9). We note that the Born (sE and the soft-virtual (S^E contributions 
have been manipulated similarly to what was done for the real-emission contribution 
in eq. (3.5); although strictly speaking this is not necessary, since these terms are 
finite and therefore not involved in any subtraction, it helps to improve the numeri¬ 
cal evaluations of the integrals I E and d§° UT ^. On the other hand, the remainders 
of the initial-state collinear subtraction (S^E only appear in /§ IN \ since /g OUT ^ is 
associated with final-state emissions. We have also introduced two two-body phase- 
space parametrizations db>E an d d0EE which are analogous to their three-body 
counterparts. Finally, we have introduced the notation 


= E 


(L=±) 

M 


+ s 


(i=/l,2) 

M 


(3.11) 
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which is the analogue of eq. (3.5) for MC subtraction terms. The first and second 
terms on the r.h.s. of eq. (3.11) receive contributions from eq. (2.39) and (2.40) 
respectively (i.e. from initial-state and final-state branchings). Taking into account 
the properties of the MC subtraction terms (see app. B), this implies that eqs. (3.6)- 
(3.9) are finite; in fact, final-state singularities of real-emission matrix elements and 
their corresponding counterterms are removed in eqs. (3.6) and (3.7) by iS (IN \ while 
initial-state singularities are removed in eqs. (3.8) and (3.9) by <S(° UT ). Thus, the 
same procedure as in sect. 4.5 of I can be used in order to generate the hard events 
that are given to the parton shower as initial conditions. 

As a concluding remark, we point out that the reason why the subtraction for¬ 
malism of refs. [17, 18] appears to be particularly well suited for MC@NLO imple¬ 
mentations can be read from eqs. (3.6)-(3.9). The partition of the phase space into 
collinear-like singular regions gives the FKS parton the same role as the softest par- 
ton emitted by an MC in the first branching after the generation of the hard process. 
Since as explained in I and II the first branching is the only one that matters for 
matching the MC with an NLO computation, the FKS parton and the softest parton 
emerging from the first branching in the shower are naturally paired in the definition 
of MC@NLO. Apart from guaranteeing the local cancellation of IR singularities, such 
pairing also allows a good control on the numerical stability of the result. It is also 
important to recall that, in each of the IR singular regions defined by the FKS parti¬ 
tion, there are no unnecessary NLO subtractions: the only counterterm contributing 
to the result is that relevant to the real matrix element singularity present in that 
given region. This fact is very beneficial in reducing the number of negative-weight 
events. 

4. Results 

In this section we present sample results for single-f production at the Tevatron with 
\fS = 1.96 TeV. We limit ourselves here to comparing MC@NLO predictions with 
those obtained with HERWIG and with an NLO code we have written according to the 
subtraction method of refs. [17, 18], as discussed in sect. 2.1.2. As a preliminary step, 
we have checked that our NLO results (with = /i ¥ = m t ) for the total rate and 
various t and t distributions are in excellent agreement with those of MCFM [9]. All 
of the predictions given in this section have been obtained by using the MRST2002 
default PDF set [19], and by setting rn t = 178 GeV, which result in total rates equal 
to 1.045 pb and 0.406 pb for t- and .s-channel respectively. We have rescaled HERWIG 
results to the NLO cross section, since we are only interested in the comparison of 
shapes in the case of standard MC’s. Also, we have only considered here HERWIG 
results for the /-channel contribution; studies that also involve the .s-charmel will 
be shown in a forthcoming paper. All the MC@NLO and HERWIG results (but 
not, of course, the NLO ones) include the hadronization of the partons in the final 
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state; furthermore, we forced the W emerging from the top decay to decay into a 
pair of leptons. In order to reduce as much as possible the statistical errors, we 
have generated 5 ■ 10 5 events for each MC@NLO and HERWIG run 5 . Finally, we 
stress that all of the fixed-order predictions presented here will be denoted as having 
NLO accuracy, even in the case of observables which, in the sense of perturbation 
theory, are effectively of leading order (see e.g. p^ J> below); this is consistent with 
the terminology one needs to adopt in the context of MC@NLO (sect. 2.3 of I). 



Figure 5: Comparison of MC@NLO (solid) and NLO (dotted) results. Left pane: top 
p T , for f-channel (higher peak) and s-channel (lower peak) contributions. We have checked 
that = Pr 1 • Right pane: top (left) and antitop (right) p, for f-channel (higher curves) 
and s-channel (lower curves) contributions. 


We start by considering the transverse momentum and pseudorapidity of the top 
and antitop (see fig. 5). We expect the impact of the momentum reshuffling that 
takes place during the hadronization phase in MC@NLO to be negligible on such 
observables. We also expect these observables, being sufficiently inclusive, to be 
reliably predicted by pure-NLO computations. As we see from the figure, the good 
agreement between MC@NLO and NLO confirms our expectations, and suggests that 
NNLO effects should be small. We have found that the HERWIG results are extremely 
close to the MC@NLO ones, and for this reason are not shown on the plots. As for 
all other processes previously studied, we have observed a much-improved behaviour 
from the numerical point of view when going from NLO to MC@NLO predictions, 
which is due to the fact that in MC@NLO all cancellations between large numbers 
occur at the level of short-distance cross sections, rather than in histograms as in the 
case of NLO computations. It is reassuring to see that this property holds true also 

5 Clearly, we are not suggesting to collect an integrated luminosity of 0(1) ab” 1 at the Tevatron. 
Here, we simply aim to expose the features of the two MC simulations with some precision. 
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for single-f production, which is the most involved process treated so far because of 
the simultaneous presence of initial- and final-state collinear singularities. 

We now discuss the properties of a few jet observables. For the sake of clarity, we 
limit ourselves in this discussion to considering ^-channel top events. We reconstruct 
the jets by means of the /c T -clustering algorithm [20], with d cut = 100 GeV 2 . We 
include in the clustering procedure all final-state stable hadrons 6 and photons. After 
the jets are reconstructed, we throw away the one that contains the 6-flavoured 
hadron whose parent parton is the b quark emerging from top decay, and order the 
remaining ones in transverse energy, i.e. the hardest jet is the one with the largest 
E t . 



Figure 6: MC@NLO (solid), HERWIG (dashed), and NLO (dotted) results, for the p T 
of the hardest jet (left pane), and the p T relative to the axis of the hardest jet of those 
hadrons or partons in that jet (right pane). 


We recall that we do not let the top decay in our pure-NLO computation. Also, 
we expect that some of the partons resulting from the radiation by the b quark 
emerging from the top decay in MC@NLO and HERWIG will hadronize into hadrons 
that are not clustered into the 6-jet which we throw away. Furthermore, some extra 
radiation will occur from the top line due to showering, which is not included in the 
NLO computation. Finally, those jets obtained with MC@NLO and HERWIG are at 
the hadron level, while those obtained with the NLO computation are at the parton 
level. 

In spite of these differences, there is a good agreement between MC@NLO, HER¬ 
WIG, and NLO for the p T of the hardest jet, shown in the left pane of fig. 6. This 
observable is sufficiently inclusive for this to happen, and the small differences be¬ 
tween MC@NLO and NLO at small p T are mainly due to the hadronization phase. 
On the other hand, the internal structure of the jet is very different in MC@NLO 

6 For the sake of simplicity, we force 7r°’s and all lowest-lying 6-flavoured states to be stable in 

HERWIG. 
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and HERWIG from that resulting from the NLO computation. In the right pane of 
fig. 6 we present the transverse momentum, relative to the axis of the jet, of all of 
the hadrons or partons clustered into the jet itself. At the NLO, the jet often coin¬ 
cides with a single parton, hence the sharp peak at p[ h 'l = 0. Such a peak is much 
less pronounced in the case of the MC’s, since in those cases the jet almost never 
coincides with a single hadron. On the other hand, at large p^h the MC results are 
smaller than the NLO one: this must be so, since in the final states obtained with 
MC simulations it is likely that a large-p^ hadron will be clustered into another jet. 
This is much less probable at the NLO, simply because the number of jets there is 
limited to two. It is also interesting to observe that, although very small, the effect 
of the hard emissions due to the NLO real matrix elements is visible in the tail of 
the p^h distribution, the MC@NLO result being slightly harder than the HERWIG 
one. 



10° 

10 _1 
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Figure 7: As in fig. 6, for the p T of the two-hardest-jet pair (left pane), and for the number 
of jets (right pane). 


The differences between the topologies of the final states emerging from NLO 
computations and MCRNLO and HERWIG simulations are clearly visible when we 
consider observables less inclusive than the p T of the hardest jet. In the left pane of 
fig. 7 we plot the p T of the pair of the two hardest jets. As is clear from the fact 
that MC@NLO and HERWIG have very similar shapes, which are different from the 
NLO one, the real matrix elements play a minor role here compared to the multiple 
emissions of the shower. The effects of the real matrix elements are more clearly 
visible in the tail of the distribution in the number of final-state jets (right pane of 
fig. 7), with MC@NLO predicting more events with more than two jets compared to 
HERWIG. 

It is also interesting to observe that shower effects dominate over matrix element 
ones for top-hardest jet correlations, two of which we present in fig. 8. We stress 
again here that we did not make any systematic attempt to exclude from the jet 
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Figure 8: As in fig. 6 , for the azimuthal difference between (left pane), and the p T of 
(right pane) the top-hardest jet pair. 


clustering the partons radiated by the top and its decay products, which would allow 
a closer matching between MC’s and NLO results for these correlations. This is very 
clearly shown by the left pane of fig. 8 , which presents the difference in azimuth be¬ 
tween the top and the hardest jet. While the NLO prediction is zero for Acp^ < 7 t /2 
for kinematics reasons (there is nothing in this region for the top-hardest jet pair to 
recoil against), MC@NLO and HERWIG feature a long tail which extends down to 
A(J)( t3> = 0. This is in part due to the fact that the top tends to have a much larger 
longitudinal than transverse momentum component. Thus, it is relatively easy for a 
parton, radiated by the top quark shower, to change the top transverse momentum 
by a sizable amount. The A= 0 tail is mainly populated by such Iow-Pt'* events. 
In the right pane of fig. 8 we present the p T of the top-hardest jet pair. At the NLO 
level, only 2^3 processes can contribute to the region p^ % 0 , in the configu¬ 
rations in which the two final-state massless partons are not combined into a single 
jet; for this to happen, the two partons must be well separated. Clearly, such con¬ 
figurations imply the presence of a very off-shell intermediate particle, and are thus 
disfavoured by matrix elements: the p. l t>> distribution is steeply falling. In MC@NLO 
and HERWIG, 2 —> 3 configurations result from a 2 —> 2 hard process followed by a 
parton branching 7 . Since the branching is collinear in nature, the probability of get¬ 
ting two well-separated partons is even smaller than in NLO computations. However, 
the shower usually does not stop after the first branching. Furthermore, all strongly- 
interacting particles, including the top and the b emerging from the top decay, can 
radiate. This smears very effectively the final-state momenta; we have verified that, 
in the large-p^' 7 '’ region, the hardest jet may retain a fraction of the parent parton 
momentum as small as 50%. This creates an imbalance between the top and the 

7 In MC@NLO, there are also 2 —> 3 hard processes, whose matrix elements are the same as 
those of the NLO computation. 
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hardest-jet p T which results in the much harder tails in MC@NLO and HERWIG 
relative to the NLO result. It should be stressed that such an effect is magnified by 
the steepness of the p-r ]> distribution. In terms of the total number of events, this 
is still a marginal phenomenon, which gives a negligible contribution to observables 
such as the inclusive p T of the hardest jet. We conclude by observing again that the 
real matrix elements contributions are small but visible in the differences between 
MC@NLO and HERWIG in the intermediate and larg e-p^ regions. 



Figure 9: MC@NLO (solid) and HERWIG (dashed) results for the p T of the ^-flavoured 
hadrons (except those from top decay). 


The results presented so far have shown little or no difference between MC@NLO 
and HERWIG results as far as shapes are concerned. Although larger differences could 
be seen by imposing hard transverse momentum cuts, the fact remains that at the 
Tevatron the phase-space for hard radiation is fairly limited. There are, however, 
observables that are particularly sensitive to real matrix element effects, such as the 
transverse momentum of the 6-flavoured hadrons 8 , which we present in fig. 9. This is 
because in i-channel matrix elements a b quark is almost always present in the initial 
state (up to CKM-suppressed contributions). This results in a final-state b -flavoured 
hadron which, in the case of HERWIG, acquires its transverse momentum entirely 
through the backward evolution in the shower mechanism. Such a mechanism is 
also present in MC@NLO, but there are also NLO real matrix elements in which a 
b quark has a large p T , which is inherited by the resulting 6-flavoured hadron, and 

8 6-flavoured hadrons from top decay are not included in this plot. 
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which explains the difference in the large-p^ tail between MC@NLO and HERWIG 9 . 

We conclude this section by mentioning the fact that we observe no dependence 
(within the statistical accuracy of the runs we performed) of the physical results upon 
the unphysical parameters which enter the NLO subtraction formalism, such as the 
subtraction parameters introduced in refs. [17, 18], or the exponent a introduced 
in eqs. (2.19) and (2.20). This constitutes a test of the correctness of our imple¬ 
mentation, since NLO results based on subtraction techniques are by construction 
independent of these parameters. Similarly, no dependence has been found on the 
parameters a and (3 introduced in eq. (I.A.86) and eq. (I.A.87) which control the 
behaviour of the MC subtraction terms in the soft limit, if they are restricted to their 
natural ranges (a = 0( 1), (3 = 0(0.1)). This is as expected, since variation of these 
parameters gives only power-suppressed effects. On the other hand, all of the above 
parameters do affect the number of negative-weight events, and their tuning can be 
used to limit the presence of such events (whose fraction is equal to about 15% in the 
results presented here). The parameter a has only a limited impact on the number of 
negative weights (which change by about 1% for 1 < a < 4), and its choice is mainly 
due to considerations of stability of the numerical integration, with best results for 
a = 2. In general, the accuracy of the predictions obtained with values of a larger 
than 2 (slowly) decreases with increasing a. Since the limit a —> oo corresponds to 
the ©-based implementation of the subtraction formalism, this indirectly proves that 
the implementation introduced in this paper is more convenient from the numerical 
point of view. 


5. Conclusions 

In this paper we have considered single-top hadroproduction in the context of the 
MC@NLO approach. This case is, apart from its phenomenological relevance, also 
interesting from the technical point of view, since it features both initial- and final- 
state collinear singularities, and thus has a radiation pattern different from that of 
all of the processes so far included in MC@NLO. 

We have shown that this is not a difficulty of principle, since the MC@NLO 
formalism is unchanged with respect to its definition given in ref. [14], but it entails 
a more involved procedure in the generation of the hard events that are given to 
the parton shower as initial conditions. Because this procedure is not specific to 
single-top hadroproduction, and since we have now treated all possible radiation 
patterns in MC@NLO, we are now in a position to include any new process, such as 
jet production, without the need of performing further analytical computations. 

9 For technical reasons, fig. 9 has been obtained by imposing \y^\ < 3. This cut has no impact 
for p^ > 10 GeV. 
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As in previous cases, our computation is based on the universal subtraction for¬ 
malism of refs. [17, 18]. We have used single-top hadroproduction as a test case, to 
explore an implementation of the subtraction different from that of the original pa¬ 
pers. The partition of the phase space is now achieved by means of smooth functions 
of invariants, rather than with 0 functions as was done previously. This does not 
entail any change in the analytical formulae, but helps to improve the behaviour of 
the numerical computations. There is also a conceptual difference, namely that the 
infrared singularities are now disentangled by means of damping factors, rather than 
by non-overlapping regions defined by the phase-space partition. This in turn may 
lead to the possibility of implementing alternative subtraction schemes, although 
new analytical computations would be required in such a case. 

We have not explored in this paper the phenomenological implications of our 
work, since we limited ourselves to checking that all of the observables we have 
considered show the expected behaviour in regions where the NLO computations or 
the MC simulations should be most reliable. We postpone the phenomenological 
studies, as well as the implementation of the Wt mode and spin correlations, to a 
forthcoming paper. 
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A. Kinematics 

In this section, we generalize the results of sect. 4 of II by considering the case of two 
final-state partons with unequal masses. Consistently with II, we use unbarred and 
barred symbols to denote quantities relevant to 2 —> 3 and 2 —> 2 processes respec¬ 
tively (see e.g. eqs. (2.16) and (2.41) for four-momentum assignments). Although in 
singlc-f production one of the final-state partons in 2 —» 2 processes is massless, we 
shall derive our results in the most general case 

— kl = ml. (A.l) 

We start by defining the 2^2 reduced invariants as follows 

s L = 2pi -p 2 , t L = ~2pi ■ ki , u L = —2pi ■ k 2 , (A.2) 
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with L = , fi , f‘ 2 ■ These invariants are used in the computations of the Born 

cross sections which appear in the MC cross sections expanded to NLO, hence the 
dependence on the branching leg L in eq. (A.2). We also get 

-2 p 2 -k 2 = t L + A m\ 2 , -2 p 2 -ki = u L - A m\ 2 , (A.3) 


where A m\ 2 = m ‘\ ~ rr A■ As discussed in II. the 2-^2 reduced invariants are 
functions of the invariants relevant to the 2^3 kinematics. The computations 
required to determine such functions are non-trivial; we only report the results here: 


S± = S + V 1 + v 2 , 


^/l,2 

= S, 






t± = 

S± 

1 X 2 (ti - Ml) + Xi(t 2 - u 2 ) 

A mj 2 

2 

[ 

2 Sy/xl - 

- XiX 2 Viv 2 /s 2 

2 

tfl = 

1 

1-1 

(t 2 - u x \ 

p' 

- ^A ml 2 


~2 S 

\s-w 1 J 

P2_ 


tf2 = 

1 

_c 

1-1 

fti-u 2 \ 

p' 

~ ^Am? 2 


-o 

2 

\ S - w 2 J 

Pi. 


u± = 

-s± - 

~i±, 

^/l,2 

= - 

~ S l/l,2 1 



(A.4) 
(A.5) 




X- \ 

Xl X 2 ViV 2 /s 2 J 
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(A.6) 


(A.7) 


(A.8) 
(A.9) 


where 


P 

Pi 

P2 




5 



4 sm\ 

{s - w 2 y ’ 



4 sm 2 

(s — Wi ) 2 


(A.10) 
(A.ll) 
(A-12) 


and E m\ 2 = iti\ + m 2 . The 2 —> 3 invariants that appear on the r.h.s. of eqs. (A.4)- 
(A.9) are labelled as in II; their definitions are also reported here in table 4. Equa¬ 
tions (A.2)-(A.12) give sufficient information, with tables 1-3, to compute the shower 
scales to be used in eqs. (2.39) and (2.40). 

We finally summarize the formulae for HERWIG showering variables. The case 
of initial-state emissions is identical to that studied in II, the condition m\ ^ m 2 
being irrelevant here. When parton 1 branches, the showering variables z + and £ + 
are related to the invariants as in eq. (II. 4.31) and eq. (II. 4.32): 


fi = -2 Ayy+iy (a. 13) 

z + 

-! = i(l-* + )(2-<;+). (A. 14) 
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Label 

Invariant 

Relation 

s 

2pi • P2 


ti 

-2pi • k x 


t2 

— 2p 2 • k 2 


Ml 

-2pi • k 2 


U 2 

-2 p 2 ■ k\ 


V\ 

-2pi • k 3 

—s — t\ — U\ 

V 2 

-2 p 2 ■ k 3 

-s - 1-2 — U 2 

Wi 

2k± ■ k 3 

S + t 2 + Ui — m\ + m 2 

w 2 

2 k 2 ■ k 3 

S + ti + U 2 + 777 1 — 7772 

M? 2 

(h + /C 2) 2 

S + V ! + V 2 


Table 4: Notation for 2 —> 3 kinematics. 


Using eq. (2.34) we can write the solutions explicitly: 


(0 _ 2|'l 
~ ,, 




1 + 


V 2 


s(l-z { l ] ) 


(A.15) 
(A.16) 


which are identical to eq. (II.4.33) and eq. (II.4.34) except for the different definition 
of the scale l. 

The branching of parton 2 will be described in terms of the variables and 
these can be obtained from eqs. (A.13)-(A.16) by interchanging variables v\ and v 2 . 

The formulae for final-state emissions are affected by the condition rn\ ^ m 2 - 
When the parton with momentum k\ branches, eq. (II.4.23) and eq. (II.4.24) still 
formally hold 


with 


= 22/i (1 - 


C/i - (1 - Zfi) 


z fi)^E 0 , 

1 + (1 ~ Z h£.h)/&^ 

1 + Pi 


5 


Pi = \j l ~ (wi +m\)/El, 

f _ (2s - (s - Wi)£ 2 )w 2 + (s - Wi) [(wi + W 2 )p 2 ~ ARlfi] 
(s - Wi)f3 2 [2s - (s - w 1 )e 2 + (s — wi)p 2 ] 

£2 = 1-^k. 

S — Wi 


(A. 17) 
(A. 18) 


(A.19) 
(A.20) 

(A.21) 
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It is apparent that eq. (A.20) coincides with eq. (II.4.27) when mi = m 2 (i.e. £2 = 1 ). 
Solving eqs. (A. 19) and (A.20) we obtain 


z 


(0 

h 


i 


(0 

h 


1 - U h 


Wl 

2 ( 1 +m 



(A.22) 
(A.23) 


which are identical to eq. (II.4.28) and eq. (II.4.29) except for the different definition 
of the scale l. 

The branching of the parton with momentum k 2 can be treated along the same 
lines. The showering variables and will be obtained from eqs. (A.22) and (A.23) 
by formally interchanging labels 1 and 2. Note that in this way the quantity £1 ap¬ 
pears in the expression for (f 2 , and 


£1 — 1 + 


Am ? 2 

s — w 2 


(A.24) 


B. MC subtraction terms 


In this section, we construct explicitly the MC subtraction terms for single-t produc¬ 
tion, expressing them in terms of the variables used in the NLO computation. In 
order to do this, we start by introducing the phase-space parametrizations used in 
ref. [17] to deal with initial- and final-state emissions; in both cases, we integrate out 
the trivial azimuthal angles. We have 


M'”’ 


1024ir 4 

s 


P ((1 - fi)s) ti diidyid cos 9d(p , 

& 


512tt 4 2-6(1 - yj ) 


6 d&dyidyjdipj , 


where /3(s) is given in eq. (A. 10), and 10 


2(1 - ml/s - 6) 
2 - 6(1 - Vj ) 


(B.l) 

(B.2) 


(B.3) 


The variables labelled with index i refer to the FKS parton (see eq. (FKS.4.3)), and 
those labelled with index j refer to the massless final-state parton that can become 
collinear to the FKS parton (see eq. (FKS.4.57)). Note that 6 is related to the 
variable x used in I and II by the following equation 


x = l -6 • 

10 Since this section is specific to single-t production, we set m 2 = 0 here. 


(B.4) 
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This implies that eq. (B.l) coincides with eq. (II.B.22). We rewrite the real-emission 
finite contributions to the single-t cross section (eq. (FKS.4.37) and eq. (FKS.4.65)) 
as follows: 

daf nJ) 

da% utJ) 
where 

> = (< . M° UT) = 6 4?” T) • (B.7) 

As discussed in I and II. MC subtraction terms can be obtained from the MC cross 
sections expanded to NLO. Thus, following eq. (II.B.21), in order to construct them 
we must write eq. (2.39) in the same form as eq. (B.5) (after relabeling), and eq. (2.40) 
in the same form as eq. (B.6) (after relabeling). In order to do this, we note that 
the Born cross sections that appear in the MC subtraction terms have the following 
forms (in order to simplify the notation, we neglect here most of the indices) 

da = M {b \s±,i±) dcos9 in , (B.8) 

l07T 

da = M^ b \s fa ,i fa )^^-dcosd out , (B.9) 

for initial- and final-state branchings respectively. Here is the Born matrix 

element, and the angles 9 in and 0 out have been introduced in eq. (II.B.32) and 
eq. (II.B.33) respectively. As discussed in II, it is not restrictive to obtain these 
scattering angles in the zero-angle-emission limits, which leads to 


+ 


*/ 5 


1 +Vi 


*/ 5J 


((1 - S<°> d$ m> ,(B.5) 


1 -Vi 


(a - yj )m {r) ) sS> d$r> , 


(B.6) 


9in = 9 , 


9 


out 1-Jj 1 


(B.10) 


where 9 and yj are integration variables in eqs. (B.l) and (B.2) respectively. The 
first relation in eq. (B.10) coincides with eq. (II.B.35). We also note that, in the 
zero-angle-emission limits, the (trivial) azimuthal angles generated by the showers 
can be chosen to coincide with the angles p and pj of eqs. (B.l) and (B.2) for initial- 
and final-state branchings respectively. We shall therefore insert the factors dp/(2n) 
and d<pj/(2n) in eqs. (2.39) and (2.40). We rewrite eq. (2.39) as follows 


da& 


MC 


= 1 n _ 2 V dt P{z±) d(£±,z±) 

a -yt)ti 2567T 4 m,vi) 

x M ( - b \s± 1 t±) O ((z±) 2 — (3(s±) d^dyid cos 9 d<p , 


(B.ll) 


where the first factor on the r.h.s. matches the “event” part of eq. (B.5), i.e. that 
obtained by replacing the distributions with ordinary functions. Using eqs. (B.l) 
and (A.4), we get 

fi{s±)dA l dy l d: cos 9dp = -——— dcp^ . (B.12) 
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Therefore 


da (±) 

cE (i=± ) 


1 
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1 1 
+ 


1 - Hi 1 ~ Hi 


(1 - y-)e 


hs( L=± ) 


d(j) 


(in) 

’3 


#f\ (B.13) 


MC , 


deft 


(in) 


MC 


= a -t ^rH^ A1<t,(s± ' f±)e ( ( ^ ,)2 -^)' (ai4) 


which is identical (up to notational differences) to eq. (I.A.72) and eq. (I.A.73). The 
MC subtraction terms that enter eq. (3.1) are readily obtained from eq. (B.14) using 
eq. (1.4.18): 


dE 


■(L=±) 


dcj) 


(in) 


d(x u ,x 2 i) cE^-^ 


MC 


d(xi,x 2 ) 


(B.15) 


MC 


The reduced Bjorken Us Xu and x^i are given in eq. (II.4.20) or eq. (II.4.22). 

We now turn to the case of final-state emissions, and we rewrite eq. (2.40) as 
follows 

= 1 (1 - Vi)t< 

MC (1 -VjRi 2567T 4 £ fa <9(6, %) 

x 0 (i - fg) 0 ( z fl -^=1 d&dVidyjdVj • ( B - 16 ) 


d(j^°d 


E, 


VV, 


Using eq. (B.2) we get 


d&dyidyjdipj = - —hiii — l iA d0^ OUT) . 

s 6' 

Inserting this equation into eq. (B.16), and using eq. (A.5) we obtain 

=(L=/«) ■ 


(B.17) 


do^E) 

dE 




MC 6 1 - Vj 
<L=f a ) 


(1 Vj) ^ d^ I ^ 


(out) 


(B.18) 


d<fi 


(out) 

’3 


2jl 2 ft(l Vj ) Ppj d((, z, ) Mm(s> ^ 

S 00 0<* OO Vj) 


x e (1 - &) e 4> - 



(B.19) 


Note that we directly defined £ (rather than £ as in eq. (B.13)) thanks to eq. (II.4.7). 

We have checked analytically that the MC counterterms introduced above locally 
cancel the collinear divergences of the real matrix elements. As already discussed in 
I and II, this happens in the soft limit only after angular integration. We therefore 
adopt here the same solutions as in eq. (II.B.43). As in the previous cases, we 
checked that the parametric dependence introduced in this way is totally negligible, 
as we expect from power-suppressed effects. 
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